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Abstract: In this paper we address the problem of interpolating a spline developable patch bounded by a given 

spline curve and the first and the last rulings of the developable surface. In order to complete the boundary of the 
patch a second spline curve is to be given. Up to now this interpolation problem could be solved, but without the 
possibility of choosing both endpoints for the rulings. We circumvent such difficulty here by resorting to degree 
elevation of the developable surface. This is useful not only to solve this problem, but also other problems dealing 
with triangular developable patches. 
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1 Introduction 

Developable surfaces have been used extensively 
in industry for modelling sheets of steel. These sur¬ 
faces are plane patches that have been curved by iso¬ 
metric transformations, preserving lengths of curves, 
angles and areas. They mimic the properties of thin 


sists in resorting to projective dual geometry. In this 
geometry “points” are planes and “planes” are points 
and th is is usef ul t o solv e the de velopab ility con dition 
I Bodduluri and Ravanil (1993); Pottmann and Farin 


( 199511 : Hu et al. (2012)). 


One 

which 


can 


also 


construct 


surfaces 


are 


steel plates that are transformed by cutting, rolling 

stead 

(Chalfant and Maekawa ( 

1998); 

or folding, but not by stretching or application of 

Pottmann and Wallner (1999); 1 

^eopoldseder ( 

2001); 

heat, which would raise manufacturing costs. 

Peternell ( 

2004): Liu et al. 


(2011); Zeng et al. 

Their inclusion in the NURBS formalism, how- 

(2012)). 

A nice review 

may be found in 

ever, has not been easy. The condition of devel- 

Pottmann and Wallner ( 

2001 

)• 

Applications to 


opability is a non-linear differential equation which 
translates into non-linear equations for the vertices 
of the control net of the surface. 

To our knowledge the first reference to NURBS 


ship hull design may be f ound in Kilgore J 1967); 


Perez and Suarez 

fcoodi . 


(120071 1 


Perez-Arribas et al. 


at General Motors 

Vlancewicz and Frev ( 

1992 

); 

Frev and Bindschadler 

(19931)1. One approach has 

been solving the developability condition for low 

degrees (Lang and Roschel 

19921); 

Chu and Sequin 

(200l; 

Chu et al. ( 

2004). 


A large famil y of Bezi e r deve l opabl e surfaces 


was 


obtained in Aumann ( 20031 2004 1 


Another approach to developable surfaces con- 


1 Corresponding author 

©Zhejiang University and Springer-Verlag Berlin Heidelberg 2013 


defining 

affine transformations between cells of the con¬ 
trol ne t. T his resu lt has been extended to spline 
(Fernandezdambrina §2007 )J and Bezier triangu¬ 
lar (Canton and Fernandez-Jambrina (2012)) devel¬ 
opable patches. A characterisation o f Bez ier ruled 
surfaces is found in Juhasz and Roth ( 2008h . 

In this paper we make use of the latter con¬ 
structions to find solutions to interpolation prob- 
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lem s with de velo pable surfaces. For instance, in 
Fernan dez- Jambrina! (120071) . we were able to draw a 
developable surface through a given boundary curve 
and two rulings, but we could not choose both end¬ 
points for these rulings. We would like to solve such 
an issue and al so ap ply the solut ion t o new problems. 

Following Fernandez-Jambrina! (2007). we first 
review in Section[5]the main features, definitions and 
the classification of developable surfaces, whereas in 
Section [3] we deal with the formalism of B-spline 
curves. In Section [4] we review the construction of 
spline developable surfaces grounded on linear rela¬ 
tions b etween v e rtices o f the B-spl ine net, that was 


Fernandez-Jambrina (20071)- In Section[5]we 


given m 

use that construction to provide solutions to an in¬ 
terpolation problem b etween a sp li ne cur ve and two 
rulings as in Fernandez -Ja mbri na (120071 1. Finally, 
in Section [0] we use degree elevation to provide our 
new solution to the problem of interpolating a devel¬ 
opable patch between a spline curve and segments of 
the rulings at both ends. This problem could not be 
solved with just our previous results. This solution is 
extended to triangular patches in Section [7) A final 
section of conclusions is included at the end of the 
paper. 

2 Developable surfaces 

A ruled surface patch fills the space between two 
parametrised curves c(u), d{u ), 

b{u, v) = (1 — v)c(u) + vd(u), u € [a, b ], (1) 

for v € [0,1], by linking with segments, named rul¬ 
ing s, the points on both curves with the same pa¬ 
rameter u. 

In general, the tangent plane to the ruled surface 
on a ruling is different for each point on the segment. 
Developable surfaces are the subcase of ruled surfaces 
for wh ich the tange n t plane is consta nt along each 
ruling (Struik (1988); Postnikov ( 19791 )). 

Let us compute a normal vector at each point of 
a ruled surface with the derivatives of the parametri- 
sation in Eq. [0 

b u (u,v) = (1 —v)c'{u)+vd'{u), b v {u,v) = d{u)—c{u), 

(b u x b v ) (u, v ) = ((1 — v)c'{u ) + vd'{u))x ( d{u ) — c{u)) 

which is linear in the parameter v. If we calculate it 
on both ends of the rulings, 

{b u x b v ) (u, 0) = c'(u) x (d(u) - c(u )), 


{b u x b v ) {u, 1) = d'{u) x {d{u) - c{u )), 

we learn that the three vectors c'{u), d'{u ), d{u) — 
c{u) are to be coplanary in order to have a constant 
tangent plane along each ruling of the surface. 
Proposition: A ruled surface parametrised as in 
Eq. [T| is developable if and only if the vector w {u) = 
d{u) — c{u), linking the points d{u), c{u), and the 
velocities d{u), d'{u) of the curves at these points 
are coplanary for every value of u. 

3 B-spline curves 

In this section we review the formalism of B- 
spline curves and their main properties in order to fix 
the n otation, which follows closely the one in I Far in 
( 2002i ). We may define a B-spline curve c{u) of de¬ 
gree n and N pieces on an interval u n +jv-i], 

so that the J-th piece of the curve is defined on 
an interval [u 7l+ /_ 2 ,u ra +/_i]. For this we require 
an ordered list of values of the parameter u, which 
are named knots , {uo,..., U 2 n+N- 2 \- The actual 
knots defining the intervals for each piece are the 
inner knots {u n - 1 ,..., u n +N- 1 } whereas the knots 
{dio,..., u n - 2 } at the beginning of the list (usually 
taken to be equal to u n ) and {u n+ N ,..., U 2 n+N- 2 } 
at the end (usually taken to be equal to u n . fjv-i) are 
auxiliary. 

Points on B-spline curves can be computed using 


the De Boor’s algorithm, c{u) = Cq\u), consisting 
on linear interpolations between consecutive vertices. 
For a curve of just one piece: 


4 °( u ) := 


^z+n ^ 
^z+n ^z+r—1 
U Ui-\-r—l 


r— 1) 
~Ca 


cLiV), (2) 


(u) 


^z+n ^z+r—1 

for i = 0,..., n — r, r = 1,..., n. 

A useful construction, named polarisation or 
blossom of the parametrisation of the curve, consists 
of interpolating in each step with a different value m 
of the parameter u, c[v\, ..., v n \ := [ui, ..., v n \, 


Ci[v I,...,v r \ :=■ 


^■z+n Vr 
'U'i+n ^z+r—1 


-4 1} [ui,...,w r -i] 


V r Ui- 1-7—1 r— l)r Vf,\ 

+ 7-7-c i+1 > i,...,u r -i{3) 


^z+n ^z+r—1 

With this notation, u <l> = 


,u, we have 


z times 


that c{u) = c[u <n> ). Vertices are recovered from the 
polarisation as Cj = c[zq,..., Uj+ n _i]. 
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These expressions are valid for B-spline curves 
with an arbitrary number of pieces, replacing the 
interval [u n — 1 , u n ] of the first piece by the interval of 
the piece under consideration. 

We may summarise some properties of the De 
Boor algorithm and the polarisation which are rele¬ 
vant for our purposes: 

1. The velocity of the curve is 

c'( u ) = ---(c” _1) (u) - Co~ 1} (m)) 

__ n (c[a <71 ~ 1> , a n ] - c[M <n ~ 1> ,it ri _i]) 

^n—1 

2. The polarisation c[v i,..., v n ] of the spline curve 
c(zz), is multiaffine and symmetric. That is, if 
A + /j = 1 , 

c[Aui + iiv i,... ,v n \ = Ac[ai,+ fic[v i,.. .,v n ]. 

Finally, we review two operations with B-spline 
curves which we shall need later on: 

Insertion of knots: Given a B-spline curve 
of degree n with vertices {co, and knots 

{zio,..., uk}, we can split into two the piece cor¬ 
responding to the interval [zz/, ztj+i] by inserting a 
new knot u, ui < u < zz/-j_i. The new list of knots is 
then obviously { zzq ,..., uk+i}, 


4 Spline developable surfaces 


The developability condition in Proposition 1 
may bjy readily now adapted to spline curves 


(F ernandez-Jambrinal ( 20071 )). 

To start, let us consider two B-spline curves of 
degree n and one segment over a common list of knots 
{zzo,..., zz 2 ra-i}, defined on the interval [zz n _i, zt„]. 
Their respective B-spline polygons are (co,..., c n }, 
{do> - - - , dn\ . 

We may draw a simple conclusion using the De 
Boor algorithm. Using Eq. [I] and the last iteration 
of Eq. [5] it is clear that the vectors c'(zz), d'(zt), 
d(zz) — c(zz) are coplanary if and only if the four points 
Cq _ 1 ) (zz), c™ _1) (zt), do _1 ^(zz), d+^zz) are coplanary 
(see Figure [1]). 



Fig. 1 Characterisation of developable surfaces 


Ui — Uj, ^ — 0 ? . . . , I) Uj-\-\ — U, Ui — Ui — 1, 

for i = I + 2,. .., K + 1, and, since the curve has not 
changed, the blossom provides the new sequence of 
vertices {c 0 ,... ,Cl + i}, 

Cj c[ui ,..., zb-|_ n _ 1 ], i 0,..., L T 1. 


The developability condition is then equivalent 
to the possibility of writing one of the points as a 
barycentric combination of the other ones. For in¬ 
stance, 

d” -1) (zz) = no{u)d^~ l \u) -b A 0 (zz)co _1) (zz) 

+ Ai (zz)c” _1) (zz), 


Degree elevation: Formally we may express a B- 
spline curve c(zz) of degree n as a curve of degree 
n + 1. The blossom c 1 of the degree-elevated 
is related to the original one in a simple form 

(Mil), 

n+l 

^c[zq , . . . , Vi— 1 , Vi- |_i, . . . , U n _|_i] 

c+i,..., v n +\] = JS -—-(5) 

and in the list of knots {zto,..., uk } the multiplicity 
of inner knots, from zz n _i to zt n +v-i, is increased by 
one, without modifying the auxiliary knots. 


curve 

(Farm 


with coefficients Ao(zt), Ai(zt), /zo(zz) = 1 — Ao(zt) — 
Ai(zz). 

We may rewrite this combination in another 
form, separating the terms related to each curve, 
also in a barycentric fashion, 

(1 — 4(zz))cg +zz) + 4(zz)c" +zz) 

= (1 - M(u))do _1 ^(zz) + M(zz)d" _1 ^(zz), (6) 

A(u) = - . X ^ U \ . . , M(u) = w , 1 w , , 
Ao(zz) + Ai(zz) Ao(zz) + Ai(zt) 

which just excludes the case of parallel vectors 
d” 1 ^(zz) — d(p+zz), c” 1 ^(zz) — Cg - +zz), which cor¬ 
responds to a cone. In this sense we use the word 
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generic, since the following results will be valid for 
all developable surfaces, but for this type of cone. 

Using blossoms and taking into account that 
these are multiaffine Eq. m, 


soms are multiaffine, 


, . . . , ILi-\-n—l 5 ] 


= C 


^2+n A 

'U'i+li • • • 5 ^i+n— 1 > Ui~r 

^2+n ^ i 


(1 — /l(zt))cg 1 ^(m) + A(u)c{ ^(u) 

- (1 +A(u)c[u <n - 1> ,u n } 

= c [u <ri_1> , (1 - A{u))u n -i + A(u)u n \ , 


the coplanarity condition (Eq. 0 may be written in 
a more compact expression, 


A*-Ui 

'Ui+n 

^2+n ^i 

^2+n A. f -i 

^2+n 

, „ r .. ... i 

“h C • • • , Ui-\-n\ 

Ui+ n ~ Ui 

U i+ n ~ A * , A* - Ui 

— Ci H - Q+i 5 

^i+n ^2 ^2+n r^2 


c[u < "- 1> , H*(u)] = d[u <n " 1> , AT (u)], (7) 


H*(u) = (1 - H(u))u n _i + A{u)u n , 
M*(u) = (1 — M(u))u n -1 + M[u)u n . 


This expression is valid for B-spline curves with 
arbitrary number of pieces, replacing the interval 
[itn_i,ii n ] of the first piece by the interval of the 
piece under consideration. 


The higher the degree of A*(u), M*(u ), the 
larger the number of conditions imposed by Eq. 0 
Hence, we restrict now to the case with constant 
A*, M *, whi ch pr oduc es the famil ies o f developable 
surfac es in lAumannl (|2003f ): I Fernandez-.Tarnbrina 
I 2007h . Jn this case expressions on both sides of 
Eq. 0may be viewed as parametrisations of curves of 
degree n — 1 and therefore this condition is equivalent 
to the same one for their blossoms, since a blossom 
is uniquely determined by its parametrisation: 
Theorem 1 Two B-spline curves of degree n and N 
pieces with the same list of knots {uo,..., uk } define 
a developable surface on the interval [u n — i, u n +N- 1 ] 
if their blossoms are related by 


c[v\, • • •, v n -\,A*] = d[vi,...,v n -i,M*], 

for some values A*, M*. 

We may obtain relations between the B- 
spline polygons of both curves by applying the 
previous expression to lists of correlative knots, 
{iq+i,..., taking into account that blos- 


since c, = c[u u ... ,u i+n _i]. 

Corollary 1: Two B-spline curves of degree n with 
the same list of knots {uq, ■ ■ ■ ,uk} and B-spline 
polygons {co,...,Cl}, {do, ■ ■ ■ ,dz,} define a devel¬ 
opable surface if the cells of the B-spline net of the 
surface are plane and their vertices are related by 

{ui-\-n A )Ci T (A Ui)Ci- (.1 

= {u i+n -M*)d i + (M*-u i )d i+1 , (8) 

for some values A*, M* and * = 0,..., L — 1. 

This family of spline developable surfaces has 
the advantage of being defined by linear relations 
between vertices, in spite of the non-linearity of the 
condition of null gaussian curvature. 

The data for this construction are the B-spline 
polygon {co,..., cl}, the list of knots {uo,..., uk} 
and, for instance, the first plane cell of the net, given 
by either do, d\ or do and the parameters A*, M*. 

Since this construction is based on blossoms 
of curves, it is compatible with algorithms for 
B-spline curves, grounded on blossoms, such as, 
for instance, the knot insertion algorithm for 
subdivision of B-spline curves. That is, if we 
split into two pieces the interval [ttj,uj_|_i] by 
inclusion of a new knot u, so that the new 
list is {uo,..., ui, u, Tt/+ 1 ,..., uk} and we com¬ 
pute the new B-spline polygons {co,... ,Cl+ i} and 
{d 0 ,..., dL+ 1 }, these new vertices satisfy Eq. [5] 

However, this construction is not compatible 
with degree elevation of B-spline curves. The degree- 
elevated B-spline developable surface through two B- 
spline curves does not coincide with the B-spline de¬ 
velopable surface through the corresponding degree- 
elevated curves. See, for instance, in Figure [5] a de¬ 
velopable surface and the control polygons of the 
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ds 



Fig. 2 Developable B-spline surface of 4 pieces of 
degree 2 

degree-elevated boundary curves (denoted by tildes): 
the central cell of the degree-elevated surface is not 
even planar. 

We show it explictly with a simple example: 
Example 1 Find a developable surface patch of de¬ 
gree two and just one piece, bounded by two curves, 
c(u) and d(u), with polygons, 

co = (0,0,0), ci = (3,3,0), c 2 = (4,3,0); 

do = (0,0,2), <h = (2,2,3), 
and knots (0, 0,1,1}. 

From Eq. [8] applied to the first cell of the B- 
spline net, i = 0, 

(u 2 —A*)c 0 +(A*— uo)ci = (u 2 —M*)d 0 +(M* —u 0 )di, 

with n = 2, uo = 0, U 2 = 1, we get 

(1 — A*)(0,0,0) + yl*(3,3,0) = (1 — M*)(0,0,2) 

+ M*( 2,2,3), 

and hence A* = — 4/3 and M* = —2. 

We lack the vertex d 2 , but for the second cell of 
the net, 

(u 3 -A*)ci+(A*- ui)c 2 = (u3~M*)di+(M*-ui)d 2 , 

|(3,3,0) - |(4,3,0) = 3(2,2,3) - 2d 2 , 

we conclude d 2 = (13/6, 3/2,9/2). 

If we formally elevate the degree of both curves 
to three, the list of knots extends to (0, 0, 0,1,1,1} 


and the new polygons obtained with Eq. [5j 


Co 

= £[0,0,0] = 

c[0,0] = c 0 = (0 

,0,0) 

Cl 

= £[0,0,1] = 

c[0,0] + 2c[0,1] 
3 

co + 2ci 

3 


= (2,2,0) 



C 2 

= £[0,1,1] = 

2c[0,l] + c[l,l] 
3 

2ci + c 2 

3 


= (10/3,3,0) 



C3 

= £[i,i,i] = 

c[l,l] = c 2 = (4 

,3,0) 


d 0 = J[0,0,0] = d[0,0] = d 0 = (0,0,2) 

d, = J[0,0,1] = dl °’ 01 + 2rf[U| 11 = 

3 3 

= (4/3,4/3,8/3) 

* = = 

3 3 

= (37/18,11/6,7/2) 

d' 3 = d[l,l,l] = d[l,l] = d 2 = (13/6,3/2,9/2) 

correspond to a developable surface with non con¬ 
stant A*(u) = —2 — u/ 2, M*(u) = —3 — u/2 and it 
is easy to check that the four points that form the 
second cell, Ci, c 2 , di, d 2 do not lie on a plane. 

This feature, however, will be shown to be use¬ 
ful for solving interpolation problems, as it will be 
apparent in the following sections. 


di di do 



Fig. 3 Degree-elevated developable surface of one 
piece of degree 2 


5 Interpolation of B-spline developable 
surfaces 

Let us consider the following interpolation prob¬ 
lem: 

Problem 1: Given a spline curve c(u) of degree n, N 
pieces, B-spline polygon (cq, ..., Cl} and list of knots 
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{ito,.. -,u K }, U <E [a, b\, a = u n - i, b = u n+N - i, and 
two straight lines l a and If, through the endpoints 
of c(u) with respective director vectors v, w, find a 
developable surface b(u, v) such that c(u, 0) = c(u) 
and l a and lb are the first and last rulings of the 
surface, that is, l a : c(a,v), lb : c(b,v). 

The special cas e of Bezier curves of degree n was 
solved by Aumann (2003), making use of his family 
of developable surf aces. Hi s so lution is extended to 
spline curves in Fernandez -Ja mbrinai 020071) , solving 
the recursion in Eq. [8]for the B-spline net. We review 
here this construction in order to extend it to solve 
new problems in next sections. 

We focus on the general case of crossing rulings 
l a and lb, since the particular cases of parallel or in¬ 
tersecting rulings may be solved in a simpler fashion 
resorting to cyli nders and cones respec tively. 


As in Fernandez-Jambrina (2007), the last rul¬ 


ing of the developable surface can be written in terms 
of the B-spline net of the curve c(u ), the list of knots 
and the coefficients A*, M *, 


L — l 


d L - C L =H M *„ U [ +n (d 0 - Co) 


M* - m 


i=o 

A* - M* 
M* - u L -1 


[c L — a(M*)), 


m,r*\ ~ Ul- i -pq- M* - u i+n 

a ( M -11 ~rr* - c ° 

M — uq -Hj M * — ui 

i= 1 \ j=i J 


From this expression we learn that the vec- 
tors along the first and last rulings, do — cq = crv, 
dL — cl = tw, and the vector, cl — a(M*) have to 
be linearly dependent and this will happen for any 
solution Mq of the algebraic equation 

det(a(M*) — cl, v, w) = 0. (10) 


This allows us to write the linear combination 
in terms of a basis {v, w, n}, n = v x w, 

a(Mg) = cl + crv + /3w + On, 

where the coefficients are readily obtained by 
Cramer’s rule, 


/3 = det(v, «(Mq ) - c L ,n) 
det(v, w, n) 

Since M* is fixed by the coplanarity condition 
in Eq. na if we wish, we can modify the length of the 
rulings through either a or r just with the parameter 
A *, which remains free so far, 


a 


T 


A* — Mq Vr MS-Ui 
a M*-u L -1 1 = 1 M*-u l+n 

o M$-A* 

P M*-u L - i ' 


( 11 ) 


Hence, we have solved the interpolation problem 
and we can use A* for fixing either do or dL, but we 
cannot choose both ends of the rulings. An example 
of this construction is shown in Figure [4] 



Fig. 4 Developable surface of degree 2 and 2 pieces 


The procedure for solving the problem is clear: 

1. Write the algebraic equation [TO] with the B- 
spline polygon for c{u), vectors v, w and the 
list of knots and obtain a solution Mq. For any 
value of A* the resulting developable surface will 
have c(u ) as part of the boundary and the first 
and last rulings will be straight lines with re¬ 
spective directions v, w. 

2. Fix Aq by choosing either do or c?l in Eq. [Ill 

3. Use the recursivity relation in Eq.[8]for comput¬ 
ing the vertices di for d(u). 

4. The B-spline polygons {c 0 ,...,Cl}, 
{doj ■ ■ •; <^l} form the B-spline net for the 
developable patch complying with the prescrip¬ 
tion. 


det(a(Mp) - c L ,w, n) 
det(v, w, n) 


We illustrate this with an example, which will 
be useful as a first step for following sections: 
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Example 2 Consider a spline curve of degree three 
and three pieces with B-spline polygon 

co = (0,0,0), Cl = (2,3,0), c 2 = (4,3,0), 
c 3 = (5,0,0), c 4 = (7,2,1), c 5 = (9, —1,3), 

and list of knots {0,0, 0,0.3,0.7,1,1,1}, not uni¬ 
formly spaced. For the first ruling we choose di¬ 
rection v = (0,0, 2) and for the last ruling we choose 
w = (—1,0,1). Find a developable surface patch 
bounded by c(u ) and the rulings defined by v, w. 
We calculate the determinant in Eq. [lUl 

det(a(M*) — cl, v, w) 

_ 2(M* 4 + 6.2M* 3 - 12.3M* 2 + 9.3 M* - 2.1) 

~~ M* 3 (M* -0.3)(M* -0.7) ’ 

and we ensure developability by choosing the param¬ 
eter M* as one of the real solutions of 

M* 4 + 6.2 M* 3 - 12.3M* 2 + 9.3M* - 2.1 = 0, 


which are M* = —7.91, 0.37. 

We further choose d 0 = cq + v = (0,0, 2) along 
the first ruling, which amounts to choosing er = 1 
in Eq. El to obtain the respective values of the pa¬ 
rameter A* = —6.18, 0.61. We perform the calcu¬ 
lations for the first pair of parameters, A* = —6.18, 
M* = -7.91. 

We may use now Corollary 1 to obtain the B- 
spline polygon of the other boundary curve of the 
developable patch through c(u) with prescribed rul¬ 
ings, 

7 _ (,'U'i+n A )a ( A Ui)c i+ 1 T (M Ui+n^di 

di+1 — 7TJ 

M* — Ui 

for i = 0 .. . L — 1, 


di 

d2 

d 3 


di 


d 3 


(u 3 - A*)cq + (A* - u 0 )ci + (M* - u 3 )d 0 



M* 

- Ug 




(1.56, 

2.34,2.08) 





(Ui - 

A*)a + {A* - 

ui)c 2 

+ 

(M* - 

- u 4 )di 


M* 

- Ui 




(3.09, 

2.29,2.26) 





(U5 - 

A*)c 2 + (A* - 

u 2 )c 3 

+ 

(M* - 

- u 5 )d 2 


M* 

- U 2 




(3.75, 

-0.15,2.55) 





(U6 - 

A*)c 3 + {A* - 

u 3 )a 

+ 

(. M * - 

- Ug)d 3 


M* 

- u 3 




(5.22, 

1.42,3.55) 





{u 7 - 

A*)a + {A* - 

Ui)c 5 

+ 

(M* - 

- u 7 )di 


M* - Ui 

(6.76,-1.00,5.24). 


and check that in fact d 3 lies on the last ruling since 

d 5 -c 5 = (-2.24,0.00,2.24), 

which is a vector proportional to w. The resulting 
patch is shown in Figure [5] 



Another way to look at this developable sur¬ 
face would be to split the spline curve into three cu¬ 
bic Bezier curves, {Co, Ci, C 2 , C 3 }, {C 3 , C 4 , C 5 , C 6 }, 
{Cq, Ci, Cg, Cg}, by knot insertion, 

C 0 = (0,0,0), Cr = (2,3,0), C 2 = (2.86,3,0), 

C 3 = (3.48,2.61,0), Ci = (4.3,2.1,0), 

C 5 = (4.7,0.9,0), Cq = (5.52,1.04,0.33), 

C 7 = (6.14,1.14,0.57), C 8 = (7,2,1), 

C 9 = (9,-1,3). 

If we also split by knot insertion 
the other boundary curve in three cubic 
pieces, {An D 3 , H 2 , D 3 }, {D 3 , D 3 , D 5 , D e }, 

{Dq, D 7 , D s , Dg}, by knot insertion, 

A) = (0,0,2), Di = (1.56,2.34,2.08), 

D 2 = (2.21,2.32,2.15), 

D 3 = (2.67,1.99,2.24), A 4 = (3.29,1.56,2.35), 

Dq = (3.55,0.58,2.46), Dq = (4.15,0.68,2.84), 

Dq = (4.15,0.68,2.84), D 7 = (4.59,0.75,3.12), 

D 8 = (5.22,1.42,3.55), Dg = (6.76,-1.00,5.24), 

it is easy to check that the three pieces of the com¬ 
posite ruled surface are in fact independent devel¬ 
opable surfaces on their respective intervals [0,0.3], 
[0.3,0.7], [0.7,1], with the same parameters A* = 
—6.18, M* = 0.61. The boundary rulings of these 
Bezier developable surfaces have been marked in Fig¬ 
ure [5] 
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6 Degree elevation of developable sur¬ 
faces 


We have seen how to interpolate a spline devel¬ 
opable surface bounded by a spline curve and two 
rulings, but we cannot choose both endpoints for 
such rulings. _This is a limita tion of the procedure 


Fernandez- J amb rina (120071) described in the pre¬ 


vious sections. A way to deal with this problem is to 
try to find a solution of hig her degree . 


As it is pointed out in lAumannl (120041 1. degree 
elevation may be used for enlarging a developable 
patch by modifying the length of the ruling segments 
of the patch. The idea is simple. We may modify the 
length of the director vector 

w(u) = d(u) — c(u), 

of each ruling by multiplication by a function f(u), 
w (u) = f(u)w(u ) = d(u) — c(u), 


and as a consequence the boundary of the surface 
patch changes. For instance the new second curve 
d(u) starts at do = Co + /(w n _i)(do — cq) and ends 
at d L = c L + f(u n+N -i)(d L - c L ). 

It is clear that this transformation just changes 
the patch of the developable surface that is covered 
by the parametrisation and it allows us to change 
the endpoints do and d^ of the first and last rulings. 
The only problem is that the curve d(u ) is no longer a 
spline of degree n. The simplest choice for the factor 
is an affine function f(u) = au + b, and in this case 
the new surface patch 

b(u, v) = (1 — u)c(u) + vd(u ) 

will be of degree (n + 1,1). An example is shown in 
Figure [Gl 

The next step will be the calculation of the 13- 
spline polygon of the new boundary of the extended 
surface patch. 

First, we obtain the blossom of the new 
parametrised curve, 

d(u) = (1 - f(u)) c{u) + f(u)d{u). 

The blossom is a (n + l)-afhne symmetric form 
d[uo, ..., u n ] for which 

d(u) =d[u <n+1> ]. 

Since f{u) is an affine function, it is already its 
own blossom, f[u\ = f{u). For the product h(u) = 



Fig. 6 Developable surface of degree 2 and 2 pieces 
stretched to a patch of degree 3 


f(u)d(u) it is simple to produce an (n-fl)-affine form 
h satisfying h[zi <n+1> ] = h(u), 

h[u 0 , ...,u n ] = f(u 0 )d[ui, ...,u n ], 


but this form is clearly non-symmetric. 

However, we may obtain a symmetric form just 
by permuting the argument of the function /, 


h[u 0 ,... ,u n \ 


^/(ui)d[uo,..., Ui- 1 , u i+ 1 ,..., u n \ 

i =0_ 

n + 1 


This form h is (n + l)-affine, symmetric and 
clearly h[u <n+1> ] = h(u). Hence, it is the blossom 
of the parametrisation h{u). 

We may use this result to conclude that the blos¬ 
som of d(u) is given by 


d[u 0 , • • • ,u n ] 


n 

^ ~2f( u i)d[uo , . . 
2=0 


71 + 1 


n 

520--f(ui))c[uo , . . . , Ui— 1 , Ui- 1_1 5 ... 5 Ufi\ 


2=0 


n + 1 


( 12 ) 


The degree of the curve c(u) must be formally 
elevated to n + 1 in order to complete the B-spline 
net of the surface patch of degree (n+ 1,1). It can be 
computed by taking / = 1 in the previous formula 
for d. The degree-elevated blossom for c(u ) is 


1 n 

c[u 0 , ■ - * , U n ] ■ —^ ' c[u 0 , • • ■ , Ui— 1 , Ui- )_i, . • - , Un \. 
n + 1 


i=0 
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The list of knots of the degree-elevated curves 


( Farm ( 120021 ) 1 is also modified by increasing by one 
the multiplicity of the inner knots u n - 1 ,... ,u n +N- 1 , 


We go back now to Example [5J 
Example 3 Consider a spline curve of degree three 
and three pieces with B-spline polygon 


{uq, • ■ • , 1 , U n _i , ■ • ■ • ^n+N — 1 , ^n+N— 1 , ■ • • , U K } ■ 


Then the new B-spline polygons of the curves 
c(u) and d(u) will be {co,..., Cz/}, (den • • ■, dt'}, 

Ci — c[zq, . . . , U,;_|_ n ], di — d[zij, . . . (13) 

for * = 0,..., L'. The list of knots has been renum¬ 
bered as {£t 0 ,... ,u k'} in order to have correlative 
indices. 

This construction is useful to solve the following 
interpolation problem: 

Problem 2: Given a spline curve c(u) of degree n, N 
pieces, B-spline polygon {cq, ... ,Cl} and list of knots 
{u 0 , ■ ■ ■,%'}, u € [at,6], a = u n - 1 , b = u n+N -i, and 
two points do, dz., find a developable surface b(u,v) 
such that c(u, 0) = c(u), c(a, 1) = do, c(6,1) = dz,. 
The procedure for solving this problem is clear: 


1. Write the algebraic equation [TO] with the B- 
spline polygon for c(u), the list of knots and 
vectors for the rulings Codg, c^d^ and obtain a 
solution Mq. 

2 . Fix A* 0 by choosing do in Eq. [TT] (<r = 1 , but 
r / 1 in general). 

3. Use the recursivity relation in Eq.[ 8 ]for comput¬ 
ing the vertices of d(u). 

4. Increase by one the multiplicity of the inner 
knots of the boundary curves. 


co = (0,0,0), Cl = (2,3,0), c 2 = (4,3,0), 

C 3 = (5,0,0), C 4 = (7,2,1), c 5 = (9,-1,3), 

and list of knots {0, 0,0, 0.3, 0.7,1,1,1}. For the first 
ruling we choose direction v = ( 0 , 0 , 2 ) and for the 
last ruling we choose w = (—1,0,1). Find a devel¬ 
opable surface patch bounded by c(it), an unknown 
curve d(zt) and the rulings defined by v, w, such that 
d( 0 ) = c 0 + v = ( 0 , 0 , 2 ), d(l) = c 5 + w = ( 8 , - 1 ,4). 

We already have obtained that the spline curve 
with B-spline polygon 

d 0 = (0,0,2), di = (1.56,2.34,2.08), 

d 2 = (3.09,2.29,2.26), d 3 = (3.75,-0.15,2.55), 

d 4 = (5.22,1.42,3.55), d 5 = (6.76,-1.00,5.24), 

and the same list of knots provides a developable 
surface patch with the required prescription except 
that ds lies on the final ruling, but it is not ( 8 , — 1 ,4). 
In fact, dg = C 5 + rw with r = 2.24. 

In order to shorten the surface patch so that 
the final vertex of the new boundary curve d(u) is 
( 8 , —1,4), we have to raise the degree of the curves 
from three to four. 

Increasing the multiplicity of the inner knots 0, 
0.3, 0.7,1, we get the new list of knots for the degree- 
elevated curves, 

{0,0,0,0,0.3,0.3,0.7,0.7,1,1,1,1}. 


5. Formally raise the degree of c{u) and compute 
the new B-spline vertices Cj with Eq. [5] 


6 . Choose f(u) so that /(a) = 1, f(b) = 1/r, 


/(«) 


b — u 
b — a 


lii-a 
r b — a 


(14) 


7. Use this function to compute the B-spline ver¬ 
tices di for the new boundary curve d(u) with 
Eq. m an d Eq. 1121 

8 . The B-spline polygons {co,..., £[/}, 
{do,..., dz,'} form the B-spline net for the 
developable patch complying with the prescrip¬ 
tion. 


We calculate first the B-spline polygon for c(u) 
as a curve of formal degree four with Eq. [5] The 
auxiliary points are computed in Appendix [A] 


co 


Cl 


C2 


C3 


c[ 0 , 0 , 0 , 0 ] = c[ 0 , 0 , 0 ] = ( 0 , 0 , 0 ) 

c[0,0,0,0.3] = e[0,0 ’°] + 3°[0,0,°.3] 

(1.5,2.25,0) 

c[0,0,0.3,0.3] = C ^°’ 0l ° 3 ^ ~^ C[ °’ ° 3 ’ ° 3] 

(2.43,3,0) 

c[0,0.3,0.3,0.7] 

c[0,0.3,0.3] + 2c[0,0.3,0.7] + c[0.3,0.3,0.7] 
4 


(3.79,2.78,0) 
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C 4 = 


C5 = 


c[0.3,0.3,0.7,0.7] 

c[0.3,0.3,0.7] + c[0.3,0.7,0.7] 


d 4 = d[0.3,0.3,0.7,0.7] = 


/(0.3)d[0.3,0.7,0.7] 


= (4.5,1.5,0) 


/(0.7)d[0.3,0.3,0.7] + (1 - /(0.3))c[0.3,0.7,0.7] 


c[0.3,0.7,0.7,1] 

c[0.3,0.7,0.7] + 2c[0.3,0.7,1] + c[0.7,0.7,1] 


= (5.21,0.51,0.14) 

= sm-rn* n c[0.7, 0.7,1] + c[0.7,1,1] 
c 6 = c]0.7,0.7,1,1] =--- 

= (6.57,1.57,0.79) 

o, = c[0.7,1,1,1] = 

= (7.5,1.25,1.5) 

c 8 = c[l, 1,1,1] = c[l, 1,1] = (9, —1,3). 

Now we have to move the curve d(u ) over the 
developable surface patch so that the new boundary 
curve d(u) goes through the endpoints of both rul¬ 
ings, shortening the director vector w(it) by a factor 
/(it) as in Ea.fldl 

f{u) = (1 -u) + 


(l-/(0.7))c[0.3,0.3,0.7] 


= (3.71,1.20,1.74) 


d 5 = J[0.3,0.7,0.7,1] = 


/(l)d[0.3,0.7,0.7] 


2/(0.7)d[0.3,0.7,1] + /(0.3)d[0.7,0.7,1] 
4 

(1 —;(l))c[0.3,0.7,0.7] 

4 

(l-/(0.7))c[0.3,0.7,1] 


(l-/(0.3))c[0.7,0.7,1] 


= (4.38,0.35,1.73) 


d 6 = d[ 0.7,0.7,1,1] = 


/(l)d[0.7,0.7,1] 


/(0.7)d[0.7,l,l] + (l-/(l))c[0.7,0.7,1] 


(l-/(0.7))c[0.7,l,l] 


= (5.68,1.30,2.14) 


2.24 


d 7 = d[0.7,1,1,1] = 


3/(l)d[0.7,1,1] + /(0.7)d[l, 1,1] 


Finally, we use Eq. [12] to compute the B-spline 
polygon of the new boundary curve of degree four 
that goes through the endpoints of both rulings, 


3(l-/(l))c[0.7,l,l] + (l-/(0.7))c[l,l,l] 
4 

= (6.56,1.05,2.70) 
d 8 =d[l, 1,1,1] = /(l)d[l, 1,1] 

+ (l-/(l))c[l,l,l] = (8,-l,4). 


do = d[0,0,0,0] = /(0)d[0,0,0] + (1 - /(0))c[0,0,0] 

= d 0 = (0,0,2) 


The degree-elevated B-spline net for the new 

surface patch, complying with the requirements of 
t /(0.3)d[0,0,0] + 3/(0)d[0,0,0.3] ,, , ’ =, 

di = d[0, 0,0, 0.3] = —- ——- w L — - - the example can be seen m Figure[7] 


(1 - /(0.3))c[0,0,0]3(1 - /(0)c[0,0,0.3] 


= (1.17,1.76,1.97) 
d 2 = d[0,0,0.3,0.3] = 


/(0.3)d[0,0,0.3] + /(0)d[0,0.3,0.3] 


(1 - /(0.3))c[0,0,0.3] + (1 - /(0)c[0,0.3,0.3] 


= (1.93,2.39,1.94) 
d 3 = d[0,0.3,0.3,0.7] = 


/(0.7)d[0,0.3,0.3] 


+ 


2/(0.3)d[0,0.3,0.7] + /(0)d[0.3,0.3,0.7] 
4 

(1 — /(0.7))c[0,0.3,0.3] 



+ 


(1 ~ /(0.3))c[0,0.3,0.7] 
2 

(1 — /(0))c[0.3,0.3,0.7] 


= (3.06,2.24,1.86) 


Fig. T Degree-elevation and restriction of the devel¬ 
opable surface patch in Figure [5] 


We could also have split the original curve c(u) 
in three cubic Bezier pieces and raise the degree of 
each of them to obtain curves of formally degree four 
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with control points, 


Co = (0,0,0), Ci = (1.5,2.25,0), C 2 = (2.43,3,0), 

C 3 = (3.01,2.90,0), C A = (3.48,2.61,0), 

C 5 = (4.09,2.23,0), C 6 = (4.5,1.5,0), 

CV = (4.910.93,0.08), C 8 = (5.52,1.04,0.33), 

C 9 = (5.99,1.12,0.51), Cio = (6.57,1.57,0.79), 
C n = (7.5,1.25,1.5), C 12 - (9,-1,3), 


and use the construction in Aumanni ( 2004 ) to extend 
each Bezier developable surface patch to comply with 
the prescription of endpoints, by multiplication by 
the same factor fiu). One reaches the same result 
as applying insertion of knots 0,0.3,0.7,1 to d(u), 


A> = (0,0,2), Di = (1.17,1.76,1.97), 

D 2 = (1.93,2.39,1.94), D 3 = (2.41,2.32,1.91), 
Di = (2.81,2.10,1.87), D 5 = (3.34,1.79,1.81), 
D 6 = (3.71,1.20,1.74), D 7 = (4.09,0.71,1.74), 
D 8 = (4.68,0.82,1.86), J % = (5.12,0.89,1.96), 
D w = (5.68,1.30,2.14), D n = (6.56,1.05,2.70), 
D V2 = { 8,-1,4). 

The boundary rulings of the quartic Bezier de¬ 
velopable surfaces have been marked in Figure 0 


7 Triangular developable surfaces 

We may pose another interpolation problem in 
which the first ruling collapses to a point, c(a) = 
d(a), 


b(u, v) = (1 — v)c(u ) + vd(u), u G [a, b]. 


We do not know the first ruling of the surface, 
but we may use previous constructions to compute a 
spline developable patch through the curve c(u ) and 
use d,L to fix the last ruling, 

b{u , v) = c{u) + vw(u), w (u) = d{u) — c{u). 

In order to collapse the first ruling to a point, 
we shorten the patch along the rulings, 

a 'll - /7 

b(u,v) = c(u)+vf(u)w(u), f(u) = - -, (15) 

b — a 

so that c(a, v) = cq for all v. 

We compute the velocity, 

A. V 

b u {u , v) = c'(u) + -- w(u) + vf(u)w\u ), 

b — a 

of the boundary curve d{u) at u = a, making use of 

Eq.Q] 

d'(a) = c n (a, 1) = c'(a) + ^ 

b — a 

ci — Co dg — Co 

— n -1— - 

u n - u n -1 b - a 

and from this expression we get the vertex do that is 
necessary for obtaining the velocity d'(a), 

do = Co + (b - a) (d!(a) - n— l - -^ , (16) 

\ 1 J 

Since we need to fix both do and d^ to obtain the 
developable patch b(u,v), the construction from the 
previous section is required and hence such a patch 
must be of degree n+1. Since c(it) is still of degree 
n, the calculation done in Eq. [12] is nonetheless valid 
whereas we keep the original vertices Co and c\. Fi¬ 
nally, shortening the surface patch as in Eg. 1151 with 
/(it) produces a triangular patch of degree n + 2. 

Summarising, the solution of this problem is re¬ 
duced to the one of Problem 2: 


The resulting developable patch is triangular in the 
sense that it is bounded by two curves and just one 
straight segment. Instead of the first point of the un¬ 
known curve of the boundary, we may give as datum 
its initial velocity d'(a). 

Problem 3: Given a spline curve c(u) of degree n, N 
pieces, B-spline polygon {co,.. •, Cl} and list of knots 
f7to, • - • , U G [d, 5], CL — T£ n _l, b — 7A n _(_jV—1, 
a point di, and a vector d'(a), find a triangular 
developable surface b(u,v) through c(it), such that 
c(ti,0) = c{u), c(a,v) = co for all v, c(b, 1) = d^, 
c u (a, 1) = d'(a). 


1. Calculate the vertex do and v = do — Co using 

Eq.m 

2. Write the algebraic equation [TO] with the B- 
spline polygon for c(u ), th e list ; of knots and 
vectors for the rulings codo, cz,di and obtain a 
solution Mg. 

3. Fix Ag by choosing do in Eq. [TT] (cr = 1, but 
r ^ 1 in general). 

4. Use the recursivity relation in Eq.[5]for comput¬ 
ing the vertices of d(it). 
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5. Increase by one the multiplicity of the inner 
knots of the boundary curves. 


so that developability is granted by choosing param¬ 
eter M* as a real solution of 


6. Formally raise the degree of c(u ) and compute 
the new B-spline vertices c, with Eq. [5] 

7. Choose /(it) so that /(a) = 1, f(b) = 1/r, 


8. Use this function to compute the B-spline ver¬ 
tices di for the new boundary curve d(u) with 
Eq. m and Eq. 1 1 21 


8 M* 4 + 2.6 M* 3 - 16M* 2 + 14.5M* - 3.5 = 0, 

that is M* = —1.92, 0.38. The other two solutions 
are complex. 

For having do = (0,0.5,2) on the first ruling, 
we need to take a = 1 in Eq. [TIJ The respective 
values of parameter A* are —1.16, 0.59. We choose 
the first pair of parameters for our calculations, Aq = 
— 1.16, Mq = 0.59. We calculate next the B-spline 


9. Increase by one the multiplicity of the inner 
knots of the boundary curves. 

10. Formally raise the degree of c(u) and compute 
the new B-spline vertices Ci with Eq. [5] 

11. Use a function /(u) = it to shrink the first ruling 
to a point and compute the B-spline vertices di 
for the new boundary curve d(u) with Eq. [HI] 
and Eq. 1121 

12. The B-spline polygons {co,..., Cl'}> 
{do, • ■ •, du} form the B-spline net for the 
triangular developable patch complying with 
the prescription. 

Example 4 Consider a spline curve of degree three 
and three pieces with B-spline polygon 

co = (0,0,0), cr = (2,3,0), ca = (4,3,0), 

c 3 = (5,0,0), C 4 = (7,2,1), C 5 = (9, —1,3), 

and list of knots {0, 0,0, 0.3, 0.7,1,1,1}. For the last 
ruling we choose direction w = (—1,0,1). Find 
a triangular developable surface patch bounded by 
c(u), an unknown curve d(u) and the ruling defined 
by w, such that d( 0) = Co, d'(0) = (20,30.5,2), 
d( 1) = c 5 + w = ( 8 , -1,4). 

First of all, we calculate the first ruling of the 
developable surface. According to Eq. [16] we need 

V = do - co = d'(0) + A( C0 - Cl ) = (0,0.5,2), 
and we calculate the determinant in Eq. 1101 
det (a(M*) — cl,v, w) 

_ 8 M * 4 + 2.6M * 3 - 16M * 2 + 14.5 M* - 3.5 
M* 3 (M* — 0.3)(M* — 0.7) ’ 


ds 



Fig. 8 Developable surface of degree 3 and 3 pieces 


polygon for the second boundary curve according to 
Corollary 1, 


di+i — 


(lli+n A )CjT(A T (iU iq_j_ n )c?2 


M* - m 


for i = 0 ... L — 1. 


do — 
di = 

d 2 = 

d3 = 

d4 = 

ds = 


(0,0.5,2) 

(it 3 - A*)c 0 + (A* - u 0 )c i + (M* - it3)d 0 
M* - u 0 

(1.21,2.39,2.31) 

(114 — A*)c\ + {A* — m)c 2 + (M* — ui)di 
M* — iti 

(2.13,2.17,3.16) 

(u 5 ~ - 4 *)c 2 + ( A* - u 2 )c 3 + (M* - n 5 )d 2 
(1.77,-0.07,4.80) 

(it6 ~ A*)c3 + (A* - it 3 )c 4 + (M* - H 6 )d 3 
M* - u 3 

(2.07,1.22,6.97) 

(117 — A*)ci + ( A* — 114)05 + (M* — 117)^4 
M* —114 

(2.92,-1.00,9.08). 
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Hence, d§ — C5 = rw, with r = 6.08. We show 
the surface patch in Figure [8] 

Next we shorten the surface patch so that the 
new boundary curve d(u) ends up at (8, —1,4). From 
the previous example we know that we are to increase 
the multiplicity of the inner knots by one, 

{0,0,0,0,0.3,0.3,0.7,0.7,1,1,1,1}, 
and formally raise the degree of c{u) to four, 
c 0 = (0,0,0), ci = (1.5,2.25,0), c 2 = (2.43,3,0), 

c 3 = (3.79,2.78,0), c 4 = (4.5,1.5,0), 
c 5 = (5.21,0.51,0.14), c 6 = (6.57,1.57,0.79), 


c 7 = (7.5,1.25,1.5), q, = (9,-1,3), 

and shorten the director vector w(u) by a factor f(u) 
as in Eq. [HO 




d 4 = d[0.3,0.3,0.7,0.7] = 


/(0.3)d[0.3,0.7,0.7] 


/(0.7)d[0.3,0.3,0.7] + (1 - /(0.3))c[0.3,0.7,0.7] 


(l-/(0.7))c[0.3,0.3,0.7] 


= (2.97,1.26,2.37) 


J 5 = d[0.3,0.7,0.7,1] = 


/(l)d[0.3,0.7,0.7] 


+ 


+ 


2/(0.7)d[0.3,0.7,1] + /(0.3)d[0.7,0.7,1] 

4 

(1 — /(l))c[0-3,0.7,0.7] 

4 

2(1 - /(0.7))c[0.3,0.7,1] + (1 - /(0.3))c[0.7,0.7,1] 


= (3.64,0.39,2.34) 
d 6 =d[0.7,0.7,l,l] = 


/(l)d[0.7,0.7,1] + /(0.7)d[0.7,1,1] 


(1 - /(l))c[0.7,0.7,1] + (1 - /(0.7))c[0.7,1,1] 


= (5.20,1.37,2.48) 


J 7 = d[0.7,1,1,1] = 


3/(l)d[0.7,l,l] + /(0.7)d[l,l,l] 


3(1 - /(l))c[0.7,1,1] + (1 - /(0.7))c[l, 1,1] 


= (6.26,1.15,2.87) 

d 8 = d[l, 1,1,1] = /(l)d[l, 1,1] + (1 - /(l))c[l, 1,1] 

= ( 8 ,- 1 , 4 ), 


so that the new boundary curve d{u) has degree four 
and B-spline polygon using Eq. 1121 given by 

do = d[ 0 , 0 , 0,0] = /(0)d[0,0,0] + (1 - /(0))c[0,0,0] 

= d 0 = (0,0.5,2) 

d! = d[0,0,0,0.3] = /(0-3)d[0,0,0] + 3/(0)d[0,0,0.3] 

(1 - /(0.3))c[0,0,0] + 3(1 - /(0)c[0,0,0.3] 

+ 4 

= (0.91,1.89,2.11) 

d 2 = d'[0,0,0.3,0.3] = /(0-3) < / [0,0,0.3] + /(0)d[0,0.3,0.3] 
| (1 - /(0.3))c[0,0,0.3] + (1 - /(0)c[0,0.3,0.3] 


where the auxiliary points are computed with blos¬ 
soms in Appendix [B] The result of this restriction of 
the surface patch is shown in Figure [9] 



co 


Restriction of the developable surface patch in 


Finally, following Eq. 1151 we further trim the 
surface patch bounded by c(u) and d(u) to shrink 
the first ruling to the vertex Cq- 

Since we are raising the degree of the curves 
from four to five, we have to increase the multiplicity 


= (1.51,2.42,2.20) 
d,3 =J[0,0.3,0.3,0.7] = 


/(0.7)e?[0, 0.3,0.3] 


Fig. 9 
Figure [8] 


2/(0.3)d[0,0.3,0.7] + /(0)d[0.3,0.3,0.7] 


+ 


(1 - /(0.7))c[0,0.3,0.3] + 2(1 - /(0.3))c[0,0.3,0.7] 


(1 — / (0))c[0.3,0.3,0.7] 


= (2.39,2.24,2.37) 


4 
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of the inner knots by one, 

{0,0,0,0,0,0.3,0.3,0.3,0.7,0.7,0.7,1,1,1,1,1}. 

The curve c(u) becomes formally of degree five 
using Eq. [5] with B-spline polygon, 

c 0 = c[ 0 , 0 , 0 , 0 , 0 ] = c[ 0 , 0 , 0 , 0 ] = ( 0 , 0 , 0 ) 

c[0,0,0,0] + 4c[0,0,0,0.3] 


ci =c[0,0,0,0,0.3] = 
= (1.20,1.80,0.0) 


C 2 = c[0,0,0,0.3,0.3] = 


2c[0,0,0,0.3] 


3c[0,0,0.3,0.3] 


= (2.06,2.70,0.0) 


c 3 = c[0,0,0.3,0.3,0.3] = 


3c[0,0,0.3,0.3] 


2c[0,0.3,0.3,0.3] , 

+ ——- - ---- = (2.66,2.96,0.0) 


c 4 = c[0,0.3,0.3,0.3,0.7] = 


c[0,0.3,0.3,0.3] 


3c[0,0.3,0.3,0.7] + c[0.3,0.3,0.3,0.7] 


= (3.69,2.69,0.0) 
c 5 = c[0.3,0.3,0.3,0.7,0.7] = 


2c[0.3,0.3,0.3,0.7] 


3c[0.3,0.3,0.7,0.7] 


= (4.34,1.79,0.0) 


C 6 = c[0.3,0.3,0.7,0.7,0.7] = 


3c[0.3,0.3,0.7,0.7] 


2c[0.3,0.7,0.7,0.7] 


= (4.66,1.27,0.03) 


C 7 = c[0.3,0.7,0.7,0.7,1] = 


c[0.3,0.7,0.7,0.7] 


3c[0.3,0.7,0.7,1] + c[0.7,0.7,0.7,1] 


= (5.31,0.72,0.20) 
c 8 = c[0.7,0.7,0.7,1,1] = 


2c[0.7,0.7,0.7,1] 


3c[0.7,0.7,1,1] 


= (6.34,1.39,0.68) 


* irn^n^i i ii 3c[0.7,0.7,1,1] 

Cg c[0.7, 0.7, 1, 1, 1J p. 


2c[0-7, 1,1,1] 


= (6.94,1.44,1.07) 


4c[0.7, 1,1,1] + c[l, 1, 1, 1] 
Cio = c]0.7,1,1,1,1] =- - - 

= (7.80,0.80,1.80) 

c n =c[l,l,l,l,l] = c[l, 1,1,1] = (9,-1,3), 

and following Eq. 1151 we shrink the rulings with a 
factor /(it ) = it. The auxiliary points are computed 


using the multiaffinity property of blossoms in Ap¬ 
pendix [C] 

Making use of Eq. [12] we obtain the B-spline 
polygon of the final boundary curve d(u) of degree 
five, 


d 0 = d[ 0 , 0 , 0 , 0 , 0 ] = /( 0 )d[ 0 , 0 , 0 , 0 ] 

+ (l-/( 0 ))c[ 0 , 0 , 0 , 0 ] = c 0 = ( 0 , 0 , 0 ) 

0,0,0,0.3] 

5 

4/(0)d[0, o, o, 0.3] + (1 - /(0.3))c[0, 0 , 0 , 0 ] 
+ 5 

, 4(1 — /(0)c[0,0,0,0.3] 


= (1.20,1.83,0.12) 
d 2 =d[ 0,0,0,0.3,0.3] = 


2/(0.3)J[0,0,0,0.3] 


3/(0)d[0,0,0.3,0.3] + 2(1 - /(0.3))c[0,0,0,0.3] 


3(1 - /(0)c[0,0,0.3,0.3] 


= (1.99,2.66,0.25) 


d 3 = d[0,0,0.3,0.3,0.3] = 


3/(0.3)J[0,0,0.3,0.3] 


2/(0)d[0,0.3,0.3,0.3] + 3(1 - /(0.3))c[0,0,0.3,0.3] 


2(1 - /(0)c[0,0.3,0.3,0.3] 


= (2.50,2.86,0.40) 


d 4 = d[0,0.3,0.3,0.3,0.7] = 


/(0.7) J[0,0.3,0.3,0.3] 


3/(0.3)J[0,0.3,0.3,0.7] + /(0)d[0.3,0.3,0.3,0.7] 
5 

(1 - /(0.7))c[0,0.3,0.3,0.3] 


+ 


3(1 - /(0.3))c[0,0.3,0.3,0.7] 


+ (l-/(0))a[0.3.0.3.0.3.0.7] = ( 3 29,2.52,0.75) 


d 5 = d[0.3,0.3,0.3,0.7,0.7] = 
3/(0.3)d[0.3,0.3,0.7,0.7] 


2/(0.7)J[0.3,0.3,0.3,0.7] 


+ 


2(1 — /(0.7))c[0.3,0.3,0.3,0.7] 
5 

3(1 - /(0.3))c[0.3,0.3,0.7,0.7] 


+ (3.65,1.64,1.09) 
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de = d[0.3,0.3,0.7,0.7,0.7] = 
2/(0.3) J[0.3,0.7,0.7,0.7] 


3/(0.7)d[0.3,0.3,0.7,0.7] 


3(1 - /(0.7))c[0.3,0.3,0.7,0.7] 
5 

2(1 - /(0.3))c[0.3,0.7,0.7,0.7] 


+ (3.83,1-15,1.30) 
d-j = d[0.3,0.7,0.7,0.7,1] = 


/(l) J[0.3,0.7,0.7,0.7] 


3/(0.7)J[0.3,0.7,0.7,1] + /(0.3) J[0.7,0.7,0.7,1] 



A 

CO 


Fig. 10 Restriction to a triangular patch of the de¬ 
velopable surface patch in Figure [9] 


(1 — /(l))c[0.3,0.7,0.7,0.7] 


prescribed, 


+ 


3(1 — /(0.7))c[0.3,0.7,0.7,1] 


+ (l-/(0-3J)g[0-7,0-7,0-7,1] = (4 _ 25; Q _ 62) L?0) 


d 8 = g?[0.7,0.7,0.7,1,1] = 


2/(l)d[0.7,0.7,0.7,1] 


3/(0.7)d[0.7,0.7,1,1] + 2(1 - /(l))c[0.7,0.7,0.7. 


3(1 — /(0.7))c[0.7,0.7,1,1] 


= (5.18,1.24,2.15) 


d 9 = d[ 0.7,0.7,1,1,1] = 


3/(l)d[0.7,0.7,1,1] 


2/(0.7)d[0.7,1,1,1] + 3(1 - /(l))c[0.7,0.7,1,1] 


2(1 -/(0.7))c[0.7,1,1,1] 


= (5.77,1.30,2.47) 


di O = d[0.7,l,l,l,l] = 


4/(l)d[0.7,1,1,1] 


/(0.7)d[l, 1,1,1] + 4(1 - /(l))c[0.7,1,1,1] 
+ 5 

(l-/(0.7))c[l, 1,1,1] 

+ 5 

= (6.67,0.72,3.03) 

dn =d[l, 1,1,1,1] = /(l)d[l, 1,1,1] 

+ (1 - /(l))c[l, 1,1,1] = dn = (8, -1,4). 


d'{ 0) = n dl _ d ° =-^-(1.20,1.83,0.12) 
= (20.00,30.50,2.00). 


8 Conclusions 

We have made use of a procedure of degree ele¬ 
vation for obtaining spline developable surfaces from 
which we know the segments of the first and last 
rulings and one of the curves of the boundary. It 
consists of first solving the problem with free end¬ 
points of the rulings and then moving the resulting 
boundary curve along the rulings to match the end¬ 
points and increase the degree of the curves by one. 
This solution is also used to solve the problem of 
finding a triangular spline developable patch from 
which we know the last ruling, one of the curves of 
the boundary and the initial velocity of the other 
curve. 
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A Auxiliary points for c(u) 


And similarly for d(u), 


We perform here calculations of auxiliary 
points for the curve c{u) over the list of knots 
{0,0, 0,0.3,0.7,1,1,1} which are needed for Exam¬ 
ple [51 taking into account that blossoms are multi- 
affine Eq. ©: 


d[0, 0,0] = d 0 = D 0 = ( 0 , 0 , 2 ) 

d[0,0,0.3] =d 1 =D 1 = (1.56,2.34,2.08) 
n 7 _ n 3 

d[ 0,0.3,0.3] =D 2 = q 7 _ q d[ 0,0,0.3] 

+ Q,3] = 0-4^i+0-3^ 

= (2.21,2.32,2.15) 
n 7 _ n 3 

d[0.3,0.3,0.3] = D 3 = ' n d[0,0.3,0.3] 


+ ■ 


0.3-0 


0.7-0 
d[0.3,0.3,0.7] = 


c[ 0 , 0 , 0 ] = 00 = 0 , = ( 0 , 0 , 0 ) 

c[0,0,0.3] = ci = Ci = (2,3,0) 

0 7_03 

c[0,0.3,0.3] = C 2 = ———-c[0,0,0.3] 


+ 


0.7-0 
0.3-0 


0.7-0 
= (2.67,1.99,2.24) 

d[0,0.3,0.7] = d 2 = (3.09,2.29,2.26) 

1 — 03 

d[0.3,0.3,0.7] = D a = - —d[ 0 ,0.3,0.7] 


0.4142 -t- O. 3 H 4 
07 


c[0,0.7,0.3] = 


0.7-0 
= (2.86,3,0) 

c[0,0.3,0.7] = c 2 = (4,3,0) 

0 7 — 03 

c[0.3,0.3,0.3] = C 3 = —c[0,0.3,0.3] 


0.4ci + 0.3c 2 
07 


0.3-0 


1 - 0 

d[l, 0.3,0.7] = 0.7d 2 + 0.3d 3 


1-0 
= (3.29,1.56,2.35) 

1-07 

d[0.3,0.7,0.7] = D 5 =- —d[ 0,0.3,0.7] 


0.7-0 
0.3-0 


0.7-0 


c[0.3,0.3,0.7] = 


+ 0.7-0 

= (3.48,2.61,0) 

1 — 03 

c[0.3,0.3,0.7] = C 4 = — — c[0,0.3,0.7] 

+ ^f^c[l,0.3,0.7] 

= 0.7c 2 + 0.3c 3 = (4.3,2.1,0) 
1-07 

c[0.3,0.7,0.7] = C 5 = --— c[0,0.3,0.7] 


0.4C 2 + 0.3 C 4 
07 


1-0 

d[l, 0.3,0.7] = 0.3d 2 + 0.7d 3 


1-0 
0.7-0 


1 - 0 
= (4.7,0.9,0) 

c[0.3,0.7,1] = c 3 = (5,0,0) 
1 - 0.7 

c[0.7,0.7,0.7] = C$ = 


c[l, 0.3,0.7] = 0.3 c 2 + 0.7c 3 


1 - 0.3 
0.7- 0.3 


c[0.3,0.7,0.7] 


1-0 

= (3.55,0.58,2.46) 

d[0.3,0.7, l]=(i 3 = (3.75,-0.15,2.55) 

1 - 0.7 

d[0.7,0.7,0.7] = H 6 = — —d[0.3,0.7,0.7] 

, 0.7-0.3 jrn , nf7 ^ 0.3D 5 + 0.4D 7 

+ TVoir^ 0 ' 7 ' 0 ' 7 ' ^ = -- 

= (4.15,0.68,2.84) 

1-07 

d[0.7,0.7,1] = D r = j—— o?[0.3,0.7,1] 

07-03 i = 03^3 + 04^ 
1-0.3 1 J 0.7 

= (4.59,0.75,3.12) 

d[0.7, 1 , 1 ] = di = D s = (5.22,1.42,3.55) 
d[l, 1,1] = d 5 = D g = (6.76, -1.00,5.24). 


c[0.7,0.7,1] 


0.3C 5 + 0.4C 7 


1 - 0.3 
= (5.52,1.04,0.33) 

1-07 

c[0.7,0.7,1] = C 7 = T ——c[0.3,0.7,1] 
0.7-0.3 r 

+ iToir c i 1 ' 0 ' 7 ' = 
= (6.14,1.14,0.57) 
c[0.7,1,1] = c 4 = C 8 = (7,2,1) 
c[l,l,l] = c 5 = C 9 = (9,-1,3). 


0.7 


~B Auxiliary points for d(u) 


We compute here auxiliary points for the curve 
d(u) over the list of knots (0,0, 0,0.3,0.7,1,1,1} 
0.3c 3 + 0.4c 4 which are needed for Example HJ using the property 
0.7 of multiaffinity (Eq. 0 for blossoms: 

d[0,0,0] =do = (0,0.5,2) 
d[0,0,0.3]=di = (1.21,2.39,2.31) 
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d[0,0.3,0.3] = 


d[0,0,0.3] 


d[0,0.3,0.7] = 
d[0.3,0.3,0.7] = 


0.7-0.3 
0.7-0 

-£^jd[0,0.7,0.3] = 
= (1.61,2.30,2.67) 

-d 2 = (2.13,2.17,3.16) 
1-0.3 


0.4g?i -f- 0.3 c?2 
07 


1-0 

0.3-0 


d[ 0,0.3,0.7] 

d[l,0.3,0.7] = 0.7d 2 + 0.3d 3 


d[0.3,0.7,0.7] = 


1-0 
(2.02,1.50,3.65) 

1 - 0.7 

—j—^-d[0,0.3,0.7] 
0.7-0 


d[l, 0.3,0.7] = 0.3d 2 + 0.7d 3 


d[0.3,0.7,1] = 
d[0.7,0.7,1] = 


d[0.7,l,l]: 
d[l,l,l] = 


1-0 
= (1.88,0.60,4.31) 

-d 3 = (1.77,-0.07,4.80) 
1-0.7 


1-0.3 
0.7- 0.3 


d[0.3,0.7,1] 


d[l, 0.7,1] = 


0.3d 3 + 0.4 c?4 


1-0.3 
= (1.94,0.67,6.04) 

= d 4 = (2.07,1.22,6.97) 

= d 5 = (2.92,-1.00,9.08). 


0.7 


C Auxiliary points for c(u) 

Finally we calculate the auxiliary points 
which are necessary to formally raise the de¬ 
gree of the curve c(u) with list of knots 
(0,0, 0,0, 0.3, 0.3, 0.7, 0.7,1,1,1,1} from four to five 
using the property of multiaffinity (Eq. [5j) for blos¬ 
soms: 

c[0,0,0,0] = c 0 = (0,0,0) 
c[0,0,0,0.3] = ci = (1.50,2.25,0.00) 

c[0,0,0.3,0.3] = c 2 = (2.43,3.00,0.00) 

0 7 — 03 

c[0,0.3,0.3,0.3] = —c[0,0,0.3,0.3] 


+ 


0.7-0 
0.3 — 0 


0.7-0 
0.4c 2 + 0.3c 3 
07 


c[0,0.7,0.3,0.3] 


= (3.01,2.90,0.00) 


c[0,0.3,0.3,0.7] = c 3 = (3.79,2.78,0.00) 

0 7 — 03 

c[0.3,0.3,0.3,0.7] = Q ? _ Q c[0,0.3,0.3,0.7] 
0.3-0 r 

+ ——-c[0.7,0.3,0.3,0.7] 


0.7-0 
0.4c 3 + 0.3c4 


= (4.09,2.23,0.00) 


c[0.3,0.3,0.7,0.7] = c 4 = (4.50,1.50,0.00) 
1-07 

c[0.3,0.7,0.7,0.7] = --—c[0.3,0.3,0.7,0.7] 


1-0.3 
0.7-0.3 
1-0.3 
0.3c 4 + 0.4c5 
07 


c[l, 0.3,0.7,0.7] 


= (4.91,0.93,0.08) 


c[0.3,0.7,0.7,1] = c 5 = (5.21,0.51,0.14) 
1-07 

c[0.7,0.7,0.7,1] = —— c[0.3,0.7,0.7,1] 


+ 


1-0.3 
0.7-0.3 


1-0.3 
0.3c5 + 0.4c6 
07 


c[l, 0.7,0.7,1] 


= (5.99,1.12,0.51) 


c[0.7,0.7,1,1] =c 6 = (6.57,1.57,0.79) 
c[0.7,1,1,1] =c 7 = (7.50,1.25,1.50) 
c[l, 1,1,1] = C8 = (9, -1,3). 

And similarly for d(u), 

J[0, 0,0,0] = Jo = (0,0.5,2) 
d[0,0,0,0.3] = (h = (0.91,1.89,2.11) 
d[0,0,0.3,0.3] =d 2 = (1.51,2.42,2.20) 
07-03- 

d[0,0.3,0.3,0.3] = - - - d[0,0,0.3,0.3] 


+ 


0.7-0 

0.3-0 


0.7-0 
0.4d 2 — 0.3d 3 
0J 


d[0,0.7,0.3,0.3] 


= (1.89,2.35,2.28) 


d[0,0.3,0.3,0.7] =d 3 = (2.39,2.24,2.37) 
07-03- 

d[0.3,0.3,0.3,0.7] = ——— d[0,0.3,0.3,0.7] 
0.3 - 0 - 

+ tt^—- d[0.7,0.3,0.3,0.7] 


0.7-0 
0.4d 3 + O. 3 J 4 
0 J 


= (2.64,1.82,2.37) 


d[0.3,0.3,0.7,0.7] = d 4 = (2.97,1.26,2.37) 
1-07- 

d[0.3,0.7,0.7,0.7] =-—d[0.3,0.3,0.7,0.7] 

JL U. «j 

07-03- 

+ ^-—d[l, 0.3,0.7,0.7] 


1 - 0.3 
O. 3 J 4 + O. 4 J 5 
07 


= (3.35,0.77,2.35) 


d[0.3,0.7,0.7,1] = d 5 = (3.64,0.39,2.34) 
1-07- 

d[0.7,0.7,0.7,1] = j——d[0.3, 0.7,0.7,1] 
+ °: 7 ~°: 3 d[l, 0.7,0.7,1] 


1 - 0.3 
O. 3 J 5 + 0.4(1(5 


= (4.53,0.95,2.42) 


0.7 


0.7 
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d[0.7,0.7,1,1] = d 6 = (5.20,1.37,2.48) 
d[ 0.7,1,1,1 \ = d 7 = (6.26,1.15,2.87) 
d[l,l,l,l]=<fe = (8,-l,4). 






